Optimizing Observational Strategy for Future /gas Constraints 
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The Planck cluster catalog is expected to contain of order a thousand galaxy clusters, both 
newly discovered and previously known, detected through the Sunyaev-Zeldovich effect over the 
redshift range < z < 1. Follow-up X-ray observations of a dynamically relaxed sub-sample of 
newly discovered Planck clusters will improve constraints on the dark energy equation-of-state found 
through measurement of the cluster gas mass fraction (fgas)- In view of follow-up campaigns with 
XMM-Newton and Chandra, we determine the optimal redshift distribution of a cluster sample to 
most tightly constrain the dark energy equation of state. The distribution is non-trivial even for the 
standard wq-Ws, parameterization. We then determine how much the combination of expected data 
from the Planck satellite and fgas data will be able to constrain the dark energy equation-of-state. 
Our analysis employs a Markov Chain Monte Carlo method as well as a Fisher Matrix analysis. We 
find that these upcoming data will be able to improve the figure-of-merit by at least a factor two. 



PACS numbers: 26.35.-(-c, 98.80.Cq, 98.80.Ft 



INTRODUCTION 

Thanks to the detailed modeling of galaxy clusters af- 
forded by X-ray observations from XMM-Newton and 
Chandra, measurements of the cluster gas mass frac- 
tion {fgas ) have become a viable cosmological probe 
[l|-|j| (see also [a-Q)- Effectively a distance measure- 
ment, fgas determinations as a function of redshift are 
a geometric probe of dark energy and, when combined 
with cosmic microwave background (CMB) anisotropy 
data from the WMAP satellite |8|, they have produced 
valuable constraints on key dark energy parameters. Us- 
ing 42 clusters, for example, [1.] place a constraint on 
the dark energy equation-of-state parameter P — wqp of 
wq = —1 ±0.14 at the 68% confidence level. One can sub- 
stantially improve this constraint by applying the tech- 
nique to a larger number of clusters (e.g., ^ 400) and 
shed new light on dark energy [^. 

To reduce the impact of modeling uncertainties, one 
preferentially uses massive and dynamically relaxed clus- 
ters, because gravity dominates their gas physics and de- 
viations from hydrostatic equilibrium should be small. 
This, however, limits the number of available systems. 
Increasing the sample size of suitable systems requires 
the surveying of large volumes. The vast majority of 
samples used to date come from the ROSAT All-Sky sur- 
vey (RASS) and extensive compilations from the mission 
archive. 

Surveys based on the Sunyaev-Zeldovich (SZ) effect 
( Q , see also the reviews by [10, \lM ) deliver a powerful 
capability of building new cluster samples. Cluster SZ 
surveys now produce catalogs of several tens of clusters 
[13 . Il3| to over one hundred [lj|, and sample size will 
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dramatically increase in the near future as more data is 
analyzed and sky covered. Of particular interest in the 
present context, the Planck satellite 15] is constructing 
the first all-sky cluster catalog since ROSAT. Thanks to 
its large survey volume, Planck will increase the number 
of hot. X-ray luminous and massive systems available for 
study at intermediate redshifts, beyond the RASS limit, 
and out to z ~ 1. It is expected to detect of order a 
thousand galaxy clusters, both newly discovered and pre- 
viously known ^. The Early SZ (ESZ) cluster list of 
the highest significance Planck detections over the com- 
plete sky was recently published [14i] along with a set of 
Planck SZ results [TMI- 

In view of planned X-ray follow-up campaigns, we ex- 
amine the potential gains in dark energy constraints from 
fgas measurements over these new cluster samples. In 
particular, we quantify the optimal redshift distribution 
for target samples. To this end, we employ the standard 
Figure of Merit (FoM ) defined by the Dark Energy Task 
Force (DETF) [2fl|] for a redshift dependent eqiiation of 
state parametrized as w{a) = wq + Wa(l — a)|2l|-[23j. 

The FoM is proportional to the inverse of the area con- 
tained in the 95% confidence contour of the wq— Wa plane. 
The DETF showed that this area is equal to the area en- 
closed in the Wa. — Wp plane, where Wp is the value of w(a) 
at the pivot scale factor ap where its uncertainty a{w,a) 
is minimized [23] • It can be easily shown that the pivot 
scale factor is: 



cov{wo,Wa) 



a- 



'■{Wa) 



1 



and thus the error on w. 
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where cov{'Wo,Wa) is the covariance between wq and Wa- 
As there is httle correlation, at least in the Gaussian 



approximation, between w^ 



Wr, 



the area is simply 



proportional to the product of the standard deviations: 

a{wp) X cr(wa)- 

This motivates the definition; 



Since clusters form from matter initially distributed 
over large volumes, we expect their composition to re- 
flect that of the cosmic mean Q. The cluster baryon 
fraction /f,, defined as the total baryonic mass over the 
total cluster mass, should faithfully reflect the cosmic 
baryon fraction 



FoM= [a{wp) X a{w^)Y 



(3) 



fb — f„ 



J stars ^ ^^b/^h] 



In this paper we look for the optimal cluster redshift 
distribution that maximizes the FoM for a combined fu- 
ture CMB+ fgas data set. We consider as illustration an 
observational X-ray campaign with a total available ob- 
servational time of ^ 5Ms to measure fgas of newly dis- 
covered clusters or to improve the fgas determinations of 
previously observed clusters. We assume that with this 
amount of time it will be possible to obtain fgas mea- 
surements for ^ 100 clusters with an error of ^ 10% and 
with a maximum redshift of Zmax = 1- This scenario 
is motivated in particular by the expectation that the 
Planck clusters will be luminous X-ray objects [la], as 
it is also suggested by the preliminary follow-up of the 
ESZ [l3|. Such an increase in sample size will certainly 
auger an improvement on dark energy constraints. Our 
study includes the CMB t emp erature and polarization 
data expected from Planck [23 , together with the cluster 
fgas measurements. Using the Fisher matrix formalism, 
we forecast uncertainties on dark energy parameters and 
calculate the FoM given different cluster redshift distri- 
butions. We determine the optimal distribution as that 
which maximizes the FoM . 

Once the optimal redshift distribution has been iden- 
tified, we perform a Monte Carlo analysis on simulated 
data and examine the improvement in constraining dark 
energy parameters when measuring fgas ■ 

Our paper is organized as follows: in section |T] we 
briefly describe the fgas method and the dependence of 
the data on cosmological parameters. Sectionlllldescribes 
the Fisher Matrix method and its use for fgas and CMB 
data, while section IIIII presents the analysis of the opti- 
mal cluster redshift distribution. Finally, in section IIVI 
we present parameter constraints from the Monte-Carlo 
study. In section |Vl we discuss our conclusions. 



I. THE fgas METHOD 

We begin by briefly recalling the fgas method and its 
use for cosmological analyses. For greater detail, we refer 
the reader to [5., 2^ . 

The X-ray gas mass fraction fgas is defined as the ra- 
tio of the X-ray emitting gas mass to the total mass of 
a cluster. It is determined from the observation of the 
X-ray surface brightness and the spectrally-determined, 
deprojected gas temperature profile. As both spherical 
symmetry and hydrostatic equilibrium are assumed for 
the clusters considered [l|, the analysis is best limited to 
the hottest, dynamically relaxed clusters. 



where f stars is the fraction of baryons in stars and /g™ 
is the true value of the cluster gas mass fraction. This 
relation is supported by hydrodynamical simulations [27|- 
[29[. Turning the relation around, allowing for systematic 
and astrophysical effects [ll, 0] and expressing the cosmic 
baryon fraction through the physical densities uJb^m = 
^b,mh^, with h being the reduced Hubble parameter h — 
i7o/(100Km/s/Mpc), we write: 



f':T{z)^^K 






b{z) 



l + s{z) 



(4) 



where 7 quantifies potential departures from hydrostatic 
equilibrium, e.g., due to non-thermal pressure support; 
K is a normalization uncertainty related to instrumental 
calibration and certain modelling issues; b(z) = bo{l + 
ahz) represents a possible cluster depletion factor and 
its evolution, and s{z) = sq{1 + asz) is the stellar baryon 
fraction as a function of redshift (see, e.g., [IJ). 

When calculating fgas from the X-ray observations, 
one needs to adopt a reference cosmological model, which 
is usually taken to be ACDM, in order, e.g., to relate the 
observed flux to the luminosity of the bremsstrahlung 
emission, related to the density of the emitting gas. It 
can be shown that the mass of the gas estimated from 
the deprojected X-ray surface brightness has a Mgas fx 

5/2 
djl dependence on the angular diameter distance, while 

the estimated total mass has a dependence of the type 

M (X dA [3|- Thus, the fgas has a a d/ dependence on 
the angular diameter distance. If the reference cosmology 
is not the true one, the gas fraction found by adopting 
the reference cosmology, /g^g (the observed fraction), is 
related to the true gas fraction via: 



/gasV^J ' 



|rcf \ 
25007 



CTi^; 



qtruc \ 
^2500 J 



r dfiz) 



Jtruc 



(^) 



1.5 



2500^^7 



'2500 



(^) 



(5) 

Here, we suppose the gas mass fraction is measured 
within an angular radius ^2500 corresponding to a phys- 
ical radius r25oc[3- If the true cosmology is different 
from the reference one, the misestimate of ^2500 intro- 
duces an additional dependence on cosmology, due to 
the fact that /gas° smoothly varies with radius [27|, 



^ The radius r25oo is defined as that radius within which the mean 
density is 2500 times the critical density of the universe at the 
redshift of the cluster. 



fgas cx {r/r25oo)^ ■ It can be shown that 6*2500 depends on 
the Hubble parameter and on the angular diameter dis- 
tance as 6*2500 = r25oo/dAiz) oc (fl"(z)d^(z))"i [2]. The 
parameter 77 is thus taken into account along with po- 
tential systematic effects. The angular diameter distance 
is defined as dA{z) = (1 + z)^^Sk{x)- Sk depends on 
cosmological parameters as: 



iinxy/W^l fc>0 



Skix) 



fc = 



(6) 



■■smhxy/Wk\ fc < 



where k is the curvature of the universe and fl^h'^ = 
h'^—ujjn—uJde is the curvature parameter, with oj™ and w^e 
the physical energy density parameters in dark matter 
and dark energy, respectively. The comoving distance x 
is defined as: 



Xiz) 



dz' 



(7) 



where H{a) is the Hubble parameter that evolves with 
redshift according to Friedmann's equation: 



H{z) 



lOOKm/s/Mpc 



= u^il + zf + VLkh^il + zf+uJdegiz) 

(8) 

where g{z) determines the evolution of the dark energy 
density. The influence of dark energy enters herein. For 
a dark energy equation-of-state parameterized as w{z) = 
Wo + Wa(l — (1 + z)~^), we have: 

g(a) = e^o -3(^.(.)+i),fe 

= (1 + z)3(l+'"o+™»)e-3to„(l-(l+2:)-i) ^g^ 

Equations H] and [5] describe the dependence of the ob- 
served /g°s on the actual cosmological parameters, which 
we constrain by fitting to the observed redshift evolution 
of /gis . Unfortunately, there are several degeneracies be- 
tween parameters that cannot be broken using fgas only 
data; e.g., systematic effect parameters such as 7 or n, 
are almost completely degenerate with uji,. Thus, priors 
or additional complemetary data are required. 

We therefore include in our analysis priors on system- 
atic effects, and combine the simulated fgas data with 
mock CMB data of the Planck satellite. In fact, Planck 
is expected to measure the physical energy density of 
baryons and of dark matter to the subpercent level, as 
well as the angular diameter distance at redshift 1100 
with very high accuracy, thus providing information on 
a combination of Hq, wq, Wa and flk- These parameters 
are, however, completely degenerate with each other in 
CMB data. For this reason, fgas data, as well as other 
redshift-distance relation probes like Supernovae [30, l3l| , 
are crucial for dark energy studies. As an example. Fig. 
[T]shows the evolution of fgas for different combinations of 
the true values of Hq , wq and Wa ■ The combinations are 



chosen to provide the same angular diameter distance at 
the redshift of recombination, i.e., to be completely de- 
generate in CMB data. Is is evident from this figure that 
the fgas data breaks this degeneracy. 
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FIG. 1. Evolution of fgas with redshift for different values of 
the true cosmological parameters. The blue solid line shows 
the evolution for true cosmology parameters equal to the ref- 
erence values, as reported in Tab. [Jand systematics as in Tab. 
im The dashed red line is for wo — —0.975 and Wa ~ —0.09 in 
the true cosmology, while the dot-dashed is for wo ~ —0.975 
and ho = 0.71. Each of these three sets of parameters provide 
similar angular diameter distances at high redshifts. 



II. FISHER MATRIX ANALYSIS 

We expect the Planck cluster catalog to significantly 
increase the number of massive systems known at inter- 
mediate redshifts jl6| . and hence the number of suit- 
able sources for fgas studies. The context guiding our 
present study is the possibility of observing ~ 100 clus- 
ters with XMM-Newton and Chandra from the newly- 
detected Planck clusters. We therefore aim to quantify 
the gain in cosmological constraints achievable with such 
a data set. To that end, we employ the Fisher Matrix 
methodology (see e.g., 132|). We will also study con- 
straints from fgas measurements combined with the CMB 
anisotropy data expected from Planck. 

For the cosmological parameters, we adopt the physical 
baryon and CDM densities, ujb = ilbh^ and Wc = ^ch^, 
the Hubble parameter Hq and the parameters describing 
the dark energy equation-of-state (EoS). We consider a 
constant EoS, w, as well as an evolving EoS of the form 
wq and Wa,: w{z) — wo + {1 — {1 + z)~^)wa. Furhermore, 
we examine both the cases of a flat universe and universes 
with free curvature flk- 

When simultaneously fitting for fgas and CMB 
data, we add in addition the normalization, As{k = 



Parameter 


Fiducial 


Ho 


71.9 Km/s/Mpc 


Qh 


0.02273 


Clc 


0.1099 


^k 





Mo 


-1 


Ws. 





T 


0.089 


Us 


0.963 


As 


2.41 X 10"'-* 



TABLE I. Fiducial WMAP5 marginalized cosmological pa- 
rameters [3311. 



Parameter 






Calibration/Modeling 


K 


l.± 0.10 


Gaussian 


Non-thermal pressure 


7 


0.9 < 1< 1.1 


Uniform 


Gas depletion 


bo 


0.66 < 0.83 < 1 


Uniform 


biz) = bo{l + abz) 


Ch 


-0.1 < 0. < 0.1 


Uniform 


Stellar mass 


So 


(0.162 ±0.049)ft?s' 


Gaussian 


s{z) = so(l -1- asz) 


as 


0.2 < < 0.2 


Uniform 


Angular correction 


V 


0.214 ±0.022 


Gaussian 



TABLE II. Fiducial values and priors for systematic param- 
eters in fgas ■ In the Fisher Matrix analysis, these priors are 
all considered as Gaussian priors. See text for references and 
details. 



0.002 /Mpc), and scalar spectral index, rig, of the primor- 
dial fluctuations, as well as the optical depth to reioniza- 
tion, T. As fiducial model for this parameter set, we use 
the best fit from the 5-y ear release of the WMAP satel- 
lite data (WMAP5) [3|, reported in Table H We verified 
that using the best fit from the 7-year release (WMAP7) 
Q does not change any of our results. Finally, we also 
consider seven additional parameters to describe system- 
atic effects as in Eq. S] and Eq. [SJ and listed in Table 
im with the assumed priors. The priors adopted are the 
same as in [1|. They reflect the current knowledge of 
their magnitude and therefore are conservative expec- 
tations for future data. Since the Fisher analysis only 
allows Gaussian priors, for those parameters with a uni- 
form prior, we choose to interpret half of the uniform 
intervals as the 2a bounds of Gaussian priors. In sec- 
tion |IVl we perform a MCMC analysis using the correct 
priors. 

For the fgas data, the Fisher Matrix is: 



F! 



'|^"'^i^(f^[C,,]-i^%if^ (10) 



fc=i 



de,^ 



de, 



where the sum is over the number of clusters N, 



clusters 

In our case, C is 

= 10% X fgas(^fc), 



and C is the covariance matrix of fgas 

diagonal and equal to C = al, with a 

as the errors on the measured fgas at different redshifts 

are not correlated and we suppose that the errors are the 

same for all clusters. 

Dividing the redshift range into Nhins , each containing 
Ubin clusters. Equation [TU] can be expressed as: 



parameters, i.e. 



8fgaa(Zbi,.) 



is largest. In fact, degenera- 



cies between parameters invalidate this deduction, mak- 
ing the optimal redshift distribution depend on the pre- 
cision of the data and its degeneracy-breaking power, as 
well as on the influence of other (non fgas ) constraints. 
For the CMB data, the Fisher Matrix is: 



F, 



CMB 



dC, 



X 



V 



^^ de, 

XY I ' 



■[c^] 



,-idC} 



de. 



(12) 



where C^ ' denotes the CMB angular power spectra 
X,Y=TT, EE, TE, plus the effective noise power spec- 
trum of Equation [13] below, and Cf^ is the covariance 
matrix of the C^ ' as defined in [3J] . We simulate angu- 
lar power spectra using the camb code[3S] adopting the 
expected values for the 143 GHz channel of the Planck 
satellite [23 ■ The experimental noise on the angular 
power spectrum is described as in [36[: 



Ne = 



with \ 



/iK-rad 



exp 



l)(^FWHMAad)^ 
8 In 2 



, (13) 



= 3 X 10~'^iJ.K^ as the temperature noise level 
(we consider a factor v2 larger for polarization noise) 
and ^FWHM = 7' for the beam size [23 ■ We neglect the 
contribution of primordial B-mode polarization spectra 
and the lensing signal and we take fsky = 1 for the sky 
coverage. We however verified that assuming a smaller 
sky coverage, e.g. fsky ~ 0.65, lowers the FoM by less 
than < 10%. 



77'^gas 

^^3 - 



E 



nhi 



Clgas (^ '^bin ) ^ ^Igas I, ^bin } 



de. 



de, 



(11) 



This equation shows explicitly the dependence of the fgas 
Fisher Matrix on the redshift distribution of the clusters, 
i.e., the number ribin of clusters in each bin. This is what 
we want to optimize. 

It is interesting to note that the tightest contraints on 
cosmological parameters do not necessarly follow from 
clusters at redshifts where the variation of fgas with the 



III. SEARCH FOR THE OPTIMAL 
DISTRIBUTION 

We calculate the total Fisher Matrix to estimate the 
cosmological parameters by summing the CMB and fgas 
Fisher Matrices: 

F^°\nu^) = F^"^^ + F'--{nu^) 

The number of clusters in each redshift bin ribin are the 
free parameters in our analysis. We search for the redshift 



distribution nbm that provides the highest value for the 
FoM , calculating the latter from the standard deviations 
CTji — {F'^°*)ii^ for each tested redshift distribution and 
using a maximization algorithm '37']. 

We first optimize the cluster distribution for the case 
of linear evolution of the dark energy equation-of-state in 
a flat universe, i.e., when w{z) — wq + {1 — (1 + z)~^)i(;a, 
constrained using a combination of CMB and fgas data. 
In Table IIIIl we report the optimal distribution for 100 
clusters in 10 bins over the redshift range z = — 1. This 
maximizes the FoM to a value of 15.2. 



Redshift Interval 


Number of clusters nti„ 


0.0- 0.1 


38 


0.1 - 0.2 





0.2 - 0.3 





0.3 - 0.4 


45 


0.4-0.5 





0.5 - 0.6 





0.6 - 0.7 





0.7- 0.8 





0.8 - 0.9 





0.9 - 1.0 


17 



TABLE III. Optimal distribution for 100 clusters in the red- 
shift range z = — 1 for the (too — Wa,i^k =0) model. 



The distribution presents a peak in the first bin at the 
lowest redshift, a second peak at an intermediate redshift 
oi z ^ 0.35 and then a third smaller peak in the last bin. 
This distribution improves the figure of merit by 30% 
compared to a simple uniform distribution, which has a 
FoM of 11.3. 

The shape of this optimized distribution is determined 
by the various factors. The expected Planck CMB data 
provide excellent contraints on the single values of ujb and 
LOc at subpercent level. However, they can only constrain 
the combination of Hq^ wq, Wa and Vl^ (but we are here 
first considering the case of a flat universe) because of 
the geometric degeneracy. Thus, the fgas distribution 
is optimized to break at best the degeneracy between 
these parameters, as well as the additional degeneracies 
introduced by the presence of fgas systematics. 

Suppose, at first, that systematics are perfectly known. 
In this case, measuring the fgas at low redshift permits 
us to measure Hq, as the effect of wq and Wa at these 
redshifts is negligible. This can be easily shown from Eq. 
[5l which in the limit of small z becomes: 



lim fl 



\{z; 



qref 

^2500 



) = /gir(^) 



rjtrue 

£o 



1.5 



(14) 



Thus, a measurement of fgas at the smallest possible red- 
shift separates the effect of H^ from wo,Wa, provided that 
Wf, and ujc are sufficiently well measured by the CMB 
data. 

Consider now the effect of redshift-independent sys- 
tematics. They can change the value of the fgas at low 



redshifts, thus being degenerate with Hq. Measuring fgas 
at high redshifts can help disentangle these effects, given 
that the redshift-dependent systematic parameters are 
limited by priors. This explains why the optimal distri- 
bution places clusters in the highest redshift bin. 

Finally, the central peak isolates the effects of wq and 
Wa- This peak is expected to be near to or at slightly 
higher redshift than the equivalence between dark energy 
and matter in the fiducial true cosmology, i.e., around 
Zeq ~ 0.42. At this redshift, the variation of fgas on true 
values of wq and Wa differs from the reference ones is 
maximized, as also shown in Fig. [TJ 

There are two issues with this result. The first is that 
the optimal distribution in Tab. IIIIl peaks in a bin cor- 
responding to redshifts smaller than z = 0.42. This is 
due to systematic effects, as their uncertainties grow with 
redshift, so measurements at low z can better constrain 
the parameters of interest. The second is the fact that 
we obviously do not know the true value of z^q] however, 
recent data (see e.g. [3^) show that the redshift of equiv- 
alence is between 0.35 < z < 0.55 at 68% c.l., in the case 
of a dark energy equation of state as considered here. 
We verified that even if the true Zeq is different from the 
fiducial one, within the current constraints, the optimal 
distribution of Tab. IIIIl still provides a FOM that is very 
close (within ^ 5%) to the maximum value. 

Thus we conclude that the use of the optimal distribu- 
tion of Tab. mil provides a FOM very close to the maxi- 
mum obtainable, even when the fiducial true cosmology 
is changed. 

We also perform several tests to check the optimal dis- 
tribution in a wide range of cases: 

- We check the stability of the shape of the opti- 
mal distribution by testing that small changes in 
the distribution still provide high FoM values. We 
perform this test by finding the distributions that 
provide a FoM within 5% of the maximum value. 
These distributions are calculated by sampling the 
FoM on a regular grid of the nun parameters, us- 
ing a step of Aubin = 10. This check confirms that 
all the distributions corresponding to a high value 
of the FoM are similar to the one in Table IIIIl i.e. 
have most of the clusters in the first, last and in- 
termediate bin around z ~ 0.42. 

This is shown in Fig. [2l where we plot the opti- 
mal distributions for the (wq , Wq , ^fc — 0) model, 
and the distributions with FoM within 5% of the 
maximum FoM obtainable with the optimal distri- 
bution. 

- We check that the optimal distribution has roughly 
the same behavior when different redshift ranges 
are considered, i.e., when the maximum redshift 
assumes a value between Zmax — 0.9 and Zmax ~ 3. 
This check has a purely theoretical interest for 
the higher redshifts, as it seems very unlikely that 
present and near-future experiments will detect 




FIG. 2. The black thick hne shows the optimal redshift dis- 
tribution of galaxy clusters in 10 redshift bins in the range 
z=0-l, corresponding to the highest FoM =15.2 (from the 
Fisher Matrix analysis). The blue lines show the redshift dis- 
tributions that have a FoM within 5% of the maximum value. 



many massive hot clusters at z > 1.5. Neverthe- 
less, we perform, this analysis to better understand 
the shape of the optimal distributions. The optimal 
distributions appear to all have similar character- 
istics for different redshift ranges. Indeed, all of 
them have around ~ 35% of the clusters in the first 
bin, ^ 15% in the last, and almost all the remain- 
ing clusters in the intermediate bin close to z^q. A 
minor peak (~ 5% of the clusters) also appears at 



1.15 when z„ 



> 1. 



Fig. 13] shows the constraints on the dark energy pa- 
rameters, as well as the FoM , for a {wq — Wa,^k = 
0) model, as a function of the maximum redshift 



of the redshift range 



for both the 



optimal and uniform distributions. It is interesting 
to note that the FoM obtained using a uniform dis- 
tribution increases with the maximum redshift till 
Zmax = 1.3 and then starts decreasing again. The 
gain in FoM offered by pushing the observations of 
a uniform distribution to such a maximum redshift 
is, in any case, small as the FoM obtained is only 
11.5. 

We also verify that a very similar optimal distri- 
bution is obtained when instead of CMB data we 
consider priors of flbh'^ = 0.0213 ± 0.002 39] and 
Hq = 73.8 ± 2.4 [40] in combination with the fgas 
data. In this case we find that the optimal distri- 
bution is almost the same as when using the full 
analysis with CMB+ fgas ■ This is due to the fact 
that now the clusters in the low redshift bin are 
used to constrain uic instead of Hq , as also this pa- 
rameter has a strong effect on the fgas at low z. For 
the z = 0—1 range in the {wq — Wa,^k = 0) model, 
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FIG. 3. FoM and the standard deviations of wo and Wa as 
a function of the maximum redshift. We determine these 
constraints by combining simulated Planck CMB data and 
fgas measurements for 100 clusters in the redshift range — 
Zmax- We show both the results for an optimal (blue solid 
lines) and a uniform (red dashed lines) distribution of clusters. 
These results are obtained with the Fisher Matrix analysis. 



the maximum FoM achievable is only 3.9. This is 
due the absence of the high redshift measurement 
of the angular diameter distance provided by the 
CMB. In fact, while the effect of wq and Wa on 
fgas , i-G., on the ratio of the angular diameter dis- 
tances (f^^ / d*^'^'^ , decreases with redshift some time 
after the epoch of equality, their integrated effect on 
the overall true angular diameter distance at high 
redshift is conserved. This is useful to constrain 
Wa, which has a weak effect on fgas (compared for 
example to wq). Thus, in absence of such a high 
redshift information, the FoM remains very small. 



We also study how much the systematics affect the 
optimal distribution and the FoM . The upper-left 
plot of Fig. m shows the optimal distribution of 
Tab. mil obtained assuming priors on systematics 
as in Tab. [TTl as well as the optimal distribution 
obtained assuming that systematic effects are per- 
fectly known. In this last case, the last bin dis- 
appears, as the degeneracy between Hq and the 
systematics is no longer present, while the inter- 
mediate peak of clusters shifts at the bin centered 
in z = 0.55. The maximum FoM obtainable in this 
case is obviously much higher than in the case when 
systematic uncertainties are present, namely FoM 
= 51.5. This indicates that the forecasted fgas data 
is not enough to break degeneracies between sys- 
tematics and cosmological parameters. Thus, pri- 
ors heavily affect the constraints on dark energy. 



We check that the optimal distribution in Tab. IIIII 
for the {wq — Wa,^k = Q) model provides high FoM 
also when other theoretical frameworks are tested. 
We consider in total 4 models: 

1. Linear redshift evolution of the dark energy 
equation of state in a flat universe {wq — 
Wa,^k = 0). This is the model we used to 
obtain the optimal distribution in Tab IIIII 

2. Linear redshift evolution of the dark energy 
equation of state with free curvature {wq — 

3. Constant equation of state in a flat universe 
(wo,Wa = 0,r2fc = 0). 

4. Constant equation of state with free curvature 

{wo.Wa = 0,17fe). 

We find the optimal distribution for each of these 
models, as shown in Fig. 21 In the figure, we also 
show the optimal distribution one would obtain if 
all the systematics were perfectly known. 

From the figure, it can be noticed that the distri- 
bution for the {wq — Wa,^k) model is very similar 
to the one obtained in the flat case. The reason is 
that now the last redshift bin is used to determine 
curvature, rather than systematic effects. In this 
way, however, systematics are constrained only by 
the priors, lowering the possible FoM achievable, 
which is FoM = 10.2. 

Considering now the case of a constant equation 
of state in a flat universe, {wo,Wa = 0,^k — 0), 
we find the optimal distribution by minimizing the 
standard deviation of wq . In this case the optimal 
distribution presents only two peaks. Roughly half 
of the clusters are in the first bin, in order to de- 
termine the Hubble constant. The rest of the clus- 
ters are placed in an intermediate bin at z = 0.55, 
where the effect of systematics is better disantagled 
from the effect of Hq. In fact, the uncertainty on 
the systematics grows with z, due to the redshift 
dependence of the stellar mass Ug and of the de- 
pletion factor ah. However, Hq, measured at low 
redshift, is mostly degenerate with all those system- 
atic effects that are redshift independent, namely 
the combination oi K, 7, bo, and sq. Thus, it is 
convenient to measure these systematics where the 
redshift-dependent effects are still subdominant, at 
redshift lower than 1. If we now consider the case 
where systematics are perfectly known, the peak at 
intermediate redshift disappears, while some clus- 
ters are placed in the last bin. This is due to the 
fact that measurements at this redshift can now 
improve the error bars obtained from the CMB on 
LUc and LOb, and thus further decrease the degener- 
acy that these parameters have with the Hubble 
parameter at low rcdshifts in fgas data. 

Finally, we consider the case of a constant equation- 
of-state with free curvature, {wq, Wa = 0, , fife). The 



Model 


FoM 


UIQ — Way^k = 
TOO ~Wa,^k 


Uniform 

11.3 

7.7 


Optimal 

Syst. No syst. 

15.2 51.5 

10.3 18.1 


TAB.|m| 
10.3 


Model 


cr(wo) 


wo,flk = 


Uniform 

0.086 
0.11 


Optimal 

Syst. No syst. 

0.073 0.035 
0.098 0.048 


TAB.|III| 

0.079 
0.10 



TABLE IV. Constraints on FoM (or a{w) in the case of a 
constant dark energy equation of state) obtained by combin- 
ing Planck CMB data and the fgaa data for 100 clusters with 
different distributions. We show the results using a uniform 
distribution, an optimal distribution calculated in each of the 
models considered and using the optimal distribution found 
for the (mo ~ Wa,^k =0) case reported in Tab. IIIII For the 
second distribution considered, we also report the results ob- 
tainable in the case of perfectly known systematic effects ('No 
syst.'). 



optimal distribution found is similar to the one for 
the flat case. 

We then compare the FoM and (j{wo) obtained with 
the optimal distributions in the different frame- 
works, with the constraints obtained using the op- 
timal distribution of Tab. IIIII in the three other 
frameworks. The results are showed in Tab II VI It 
is evident from the table that using the optimal dis- 
tribution of Tab. IIIII also in these other cases leads 
to FoM which are very close to the ones obtainable 
by using the proper optimal distribution for each 
case. For example, for the {wo,Wa = 0,ilk = 0) 
model, the error obtained on wq using the proper 
optimal distribution is cr{wo) — 0.073, while the one 
obtained using the optimal distribution of Tab. IIIII 
yields to cr(wo) = 0.079; the variation is less than 
~ 10%. Thus we conclude that the distribution 
in Tab. IIIII can be used in a wide range of mod- 
els to obtain tight constraints on the dark energy 
equation-of-state. 



IV. MCMC ANALYSIS 

A. Method 

In order to check our previous results obtained us- 
ing Fisher matrix forecasts, we derive the constraints on 
cosmological parameters using a modified version of the 
publicly available Markov Chain Monte Carlo package 
cosmomc |4l| . We use a convergence diagnostics based on 
the Gelman and Rubin statistic performed on 4 chains, 
and the Parameterized Post-Friedmann(PPF) prescrip- 



60 


- 1 1 1 1 . 


fe 50 


- n - 

- M - 


If) 

^ 40 
o 


1 A ' ' ■ 


o 30 


1^, / \ / '; ^ 


^ 20 


A / \ ■' ' ^ 


§ 10 
z 




\\ 1 \ '\ A 



0.2 0.4 0.6 0. 
z 





6U 
50 


J 


o 


40 


■ A ^ 


O 


30 


\ /'A '" 


0) 


20 


'\ / ''' ~'\ ' ^ 




10 


A /' ~\ Ai 







\ / \ / : 



0.2 



0.4 



0.6 





bO 




tn 




3 : 


11 


50 


\t '■ 








(/) 




'\\ / \ " 


J 


40 


r\, / \ -J 


o 




: V / \ t 


"o 


30 


H / \ h 


^ 




- \ \ ' - 


_Q 


20 


^ \ / \ ' ~' 


!= 




- \ \ 1 - 


T 


10 


'- \ \ / -j 


z 





\ \ \ : \ 



0.2 



0.4 0.6 
z 



0.8 




0.2 



0.4 0.6 0. 



FIG. 4. Optimized distributions of clusters for the {wo-Wa, ^k = 0) model (top-left), for (wo-Wa, f2fc) (top-right), for (wo,Wa = 
0, Jlfc = 0) (bottom-left) and for {woWa ~ 0,^1^) (bottom-right). We show the results for the case using priors on systematics 
as in Tab [n] (black solid lines) and assuming perfectly known systematics (red dashed lines). 



tion for the dark energy perturbations [42| implemented 
by [431 in the camb code. Moreover, we use the code of 



4J] to implement the fgas data in the cosmomc code. 



We simulate datasets for the Planck experiment follow- 
ing the commonly used approach, described, for example, 
in [4^- We first generate the fiducial theoretical angular 
power spectra Cf^ , Cf^ , Cf^ with the camb code for a 
cosmological model compatible with the WMAP5 maxi- 
mum likelihood parameters reported in Table HI We then 
add the experimental noise Ni described in Eauation ll3l 

For fgas we simulate a dataset of 100 clusters over the 
redshift range z — — 1 , distributed with either the opti- 
mal configuration reported in Table Hill or with a uniform 
distribution. We assume that the true and the reference 
cosmology are the same as the fiducial model of Table HI 
A 10% error is assigned to the measurements. 

We sample a set of cosmological parameters similar to 
the one used in the Fisher Matrix analysis, adopting flat 
priors: the physical baryon and cold dark matter densi- 
ties, ojb = ^bh"^ and Wc — ^ch^, the ratio of the sound 
horizon to the angular diameter distance at decoupling, 
ds (used instead of Hq), and the curvature fifc. For the 
CMB data, we additionally consider the scalar spectral 



index ns, the overall normalization of the spectrum A at 
k — 0.002 Mpc^^, and the optical depth to reionization, 
T. Furthermore, for fgas we consider the seven system- 
atic parameters described in Section |4l with the priors in 
Table nil 

We test the four models considered in the Fisher Ma- 
trix analysis, i.e., a linear evolution of the dark energy 
equation of state wq — w^ in a flat and non-flat universe, 
and a constant value of the dark energy equation-of-state 
wq in a flat and non-flat universe. 



B. Results 

Table |Vl shows our constraints on the most interest- 
ing varied and derived parameters for all the models and 
combinations of data tested, obtained using the cosmomc 
code. Concerning the {wq — Wa,^k — 0) model, the FoM 
found with a combination of the synthetic Planck CMB 
data plus an optimal distribution of fgas data is smaller 
than the Fisher Matrix value in Table ITVl (14.4. instead 
of 15.1) due to better accounting of uniform priors and a 
more reliable treatment of the non-Gaussian probability 



distributions of the parameters. Nevertheless, the combi- 
nation of Planck CMB+fgas data with an optimal distri- 
bution still improves the FoM by ^ 26% compared to a 
uniform distribution over the same redshift range, which 
provides a FoM of 11.4. 

Moreover, the results obtained with the combination 
of synthetic Planck CMB-f optimal fgas represent a ma- 
jor improvement over current data. For instance, we 
calculated the FoM from WMAP5 CMB combined with 
the fgas data for 42 clusters measured by Allen et al.[l|. 
These datasets yield a FoM of 8.3; therefore, future data 
will improve this result by almost a factor ^ 2, as clear 
from the second line of Tab. |Vl An even more impressive 
improvement is obtained if in addition to the synthetic 
Planck CMB data and the fgas data for 100 clusters with 
an optimal distribution, we add the known 42 fgas mea- 
surements of Allen et al. In this case the FoM will in- 
crease of almost a factor 3. 

The results in Table |V] also demonstrate that the 
improvement in the determination of the dark energy 
equation-of-state will be mainly driven by more accurate 
fgas data rather than better CMB data. In fact, com- 
bining Planck CMB data with Allen et al.'s fgas mea- 
surements only provides a FoM slightly better than that 
the one obtained with WMAP5 data — nominally FoM 
=8.8. 

Fig. [5]shows the 68% and 95% likelihood contour plots 
in the wq - Wa plane, obtained using current data and 
synthetic Planck CMB plus fgas data for the 42 clusters 
reported in Allen et al. JJ plus fgas forecasted data with 
an optimal cluster distribution. Fig. [6] compares the 
constraints achievable by Planck CMB only and Planck 
CMB-f /gas forecasted data for 100 clusters with an opti- 
mal distribution. In this plot, the degeneracy-breaking 



(H 




FIG. 5. 68% and 95% likelihood contour plots on the wo 
- Wa plane for WMAP5 plus fgas data for the 42 clusters 
reported in Allen et al. [l|] (red) and synthetic Planck plus 
fgas data for the 42 clusters reported in Allen et al.[l||, plus 
fgas data for 100 clusters between redshift and 1 with an 
the optimal distribution (yellow). Here the model considered 
is (wo — Wa,i^k = 0) in a flat universe. 
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FIG. 6. Constraints for the (wo — Wa,i^k ~ 0) model from 
synthetic Planck CMB data (red) and CMB-f fgas data (blue) 
for 100 clusters with an optimal distribution between z = 
0—1. We show the marginalized one-dimensional distributions 
and two-dimensional 68% and 95% c.l. contour plots. 



power of the fgas data on the parameters that deter- 
mine the geometry of the universe, such as the dark en- 
ergy equation of state parameters wo, Wa and the Hubble 
constant, is striking. By themselves, CMB data can only 
place weak constraints on these parameters when they 
are considered simultaneously. As expected, adding fgas 
data to CMB data do not improve the constraints on 
other parameters, namely the baryon and dark matter 
densities Wb and Wc, the spectral index n^, and the opti- 
cal depth T. 

Very similar considerations can be made for a constant 
equation of state w = wq in a flat universe (51^ = 0), as 
show in Figure [T] In this case, the constraint on wq 
is (j(wq) = 0.091, which represents a 35% improvement 
compared to the results reported in Allen et al. using 
WMAP3 CMB data and 42 clusters, namely cr(wo) = 
0.14. 

Finally we repeated the analysis for also the case of a 
universe with free curvature. Fig. [8]shows the constraints 
obtained for the (wq — Wat^k) model. The main effect 
of adding free curvature is to introduce more degeneracy 
with the parameters of interest. This leads to a widening 
of the errors and a distortion of the probability distribu- 
tions. This also explains why the FoM obtained with the 
MCMC analysis in (w^-Wa, ^k) case, FoM =7.5, is 30% 
smaller than the one estimated with the Fisher Matrix, 
FoM =10.3, that assumes Gaussian likelihoods for the 
considered parameters. 



V. CONCLUSIONS 

We consider constraints on the dark energy equation- 
of-state obtained when combining the future CMB data 
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TABLE V. Figures of merit and 68% c.l. errors on cosmological parameters for different models, obtained using the cosmomc 
code. The dark energy and dark matter density parameters fide and Qm are indirectly derived. We show the results for different 
combinations of data. 'WMAP5' indicates CMB data from the 5-year release of the WMAP satellite data 33], while 'Planck' 
indicates simulated CMB data from the Planck satellite [25|]. 'fgas (opt).' indicates fgaa simualated data for 100 clusters with 
an optimal distribution between redshift and 1, while 'fgas (uni.)' is for a uniform distribution. Finally, 'Allen' indicates fgas 
real data for 42 clusters publicated in Allen et al. [l|]. 



expected from the Planck satellite and simulated fgas 
measurements for 100 clusters observed by a possible fu- 
ture X-ray campaign with X-ray satellites. We found 
the optimal cluster distribution over the redshift range 
z = — 1 to provide the tightest constraints on the dark 
energy equation-of-state. The optimal distribution im- 
proves the FoM by almost ^ 30% compared to a uniform 
distribution. We found that the optimal distribution, 
found for a linearly evolving dark energy equation of state 
in a flat universe, provides very high figures of merit also 
in other theoretical frameworks, namely free curvature or 
a constant equation-of-state. A positive characteristic of 
the optimal distribution is that clusters are placed at low 
and intermediate redshifts, while only a few should be 
observed at high redshift, where measurements require 
longer observational times. We performed our analysis 



by both using the Fisher Matrix and the MCMC meth- 
ods to forecast constraints. We found that future data 
will increase the FoM by a factor ^ 2 relative to the 
current constraints available with WMAP 5 and the fgas 
data from [1[ . Our study provides a useful guide for plan- 
ning XMM-Newton and Chandra follow-up observations 
of Planck and SPT clusters. 
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